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ABSTRACT 



The "ramp" is an often observed feature in temperature 
fluctuations during unstable atmospheric conditions. It is 
characterized by a gradual increase in temperature followed 
by a sudden drop to an ambient level. This ramp clearly 
distinguishes temperature signals from other turbulence signals 
such as velocity. Three different ramp- type atmospheric 
temperature fluctuations and their derivatives are constructed 
and statistically examined for the parameters skewness and 
coefficient of excess. These statistical values are compared 
with values obtained from actual signals. The "complex ramp 
signal" was found to represent these quite well. Skewness and 
coefficient of excess for this signal could be altered almost 
iiidependently of each other. The complex ramp is constnacted 
by inscribing nearly symmetrical consecutive triangles within 
two envelope rays of constant but different slopes such that 
the two envelope rays approximate a ramp function. 



2 



I 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

A. OBJECTIVE 8 

B. GENERAL FEATURES OF TEMPERATURE SIGNALS ... 9 

II. SIGNALS L3 

A. ATMOSPHERIC TEMPERATURE FLUCTUATIONS 13 

B. MODELS OF ATMOSPHERIC TEMPERATURE FLUCTUATIONS 13 

III. METHOD OF EXAMINATION 16 

A. SIGNAL 16 

B. DIFFERENTIATED SIGNALS 18 

C. STATISTICS APPLIED 19 

IV. RESULTS 26 

V. CONCLUSIONS 49 

VI. RECOMMENDATIONS FOR FURTHER STUDY 50 

COMPUTER PROGRAM 53 

BIBLIOGRAPHY 59 

INITIAL DISTRIBUTION LIST 61 

FORM DD 1473 62 



3 



LIST OF FIGURES 



FIGURE 1: Atmospheric Temperature and Differentiated 

Temperature Signal 

FIGURE 2; Schematic of Convective Plume According to 

Gill (1971) 12 

FIGURE 3: Growth of Ramps as Atmosphere Changes from 

Stable to Unstable 15 

FIGURE 4: Standard Ramp Model and Differentiated 

Standard Ramp Model ....o 22 

FIGURE 5; Asymmetrical Ramp Model and Differentiated 

Asymmetrical Ramp Model 23 

FIGURE 6 : Complex Ramp Model and Differentiated 

Complex Ramp Model » 0....24 

FIGURE 7: Variation of Standard Ramp Signal with 

Period 28 

FIGURE 8 : Variation of Standard Ramp Signal with 

Ramp Width 31 

FIGURE 9: Variation of Standard Ramp Signal with 

Ramp Slope 31 

FIGURE 10: Variation of Asymmetrical Ramp Signal with 

Period 33 

FIGURE 11: Variation of Asymmetrical Ramp Signal with 

Ramp Width 34 

FIGURE 12: Variation of Asymmetrical Ramp Signal with 

Back Slope 34 

FIGURE 13: Variation of Asymmetrical Ramp Signal with 

Back Slope Width 38 

FIGURE 14: Variation of Complex Ramp Signal with 

Period ...40 

FIGURE 15: Variation of Complex Ramp Signal with 

Upper Basic Slope 42 



4 



\ 



FIGURE 16: Variation of Complex Ramp Signal with 

Lower Basic Slope 44 

FIGURE 17: Variation of Complex Ramp Signal with 

Ramp Width 45 



5 



I 



LIST OF TABLES 



TABLE 


I. 




SKEWNESS AND COEFFICIENT OF EXCESS OF 
DIFFERENTIATED TEMPERATURE SIGNALS .... 


25 


TABLE 


II. 


- 


VARIATION OF DIFFERENTIATED 
WITH PERIOD 


RAMP SIGNAL 


27 


TABLE 


III. 


- 


VARIATION OF DIFFERENTIATED 
WITH RAMP SLOPE 


RAMP SIGNAL 


29 


TABLE 


IV. 


- 


VARIATION OF DIFFERENTIATED 
WITH RAMP WIDTH 


RAMP SIGNAL 


30 


TABLE 


V. 


- 


VARIATION OF DIFFERENTIATED 
RAMP SIGNAL WITH PERIOD. . . 


ASYMMETRICAL 


32 


TABLE 


VI. 


- 


VARIATION OF DIFFERENTIATED 
RAMP SIGNAL WITH BACK SLOPE 


ASYMMETRICAL 


34 


TABLE 


VII. 


- 


VARIATION OF DIFFERENTIATED 
RAMP SIGNAL WITH RAMP WIDTH 


ASYMMETRICAL 


35 


TABLE 


VIII. 


- 


VARIATION OF DIFFERENTIATED 
RAMP SIGNAL WITH BACK SLOPE 


ASYMMETRICAL 
WIDTH .... 


37 


TABLE 


IX. 


- 


VARIATION OF DIFFERENTIATED 
SIGNAL WITH PERICD .... 


COMPLEX RAMP 


39 


TABLE 


X. 


- 


VARIATION OF DIFFERENTIATED COMPLEX RAMP 
SIGNAL WITH UPPER BASIC SLOPE 


41 


TABLE 


XI. 


- 


VARIATION OF DIFFERENTIATED COMPLEX RAMP 
SIGNAL WITH LOWER BASIC SLOPE 


43 


TABLE 


XII. 


- 


VARIATION OF DIFFERENTIATED 
SIGNAL WITH RAMP WIDTH . . , 


COMPLEX RAMP 


45 


TABLE 


XIII. 


- 


VARIATION OF DIFFERENTIATED 
SIGNAL WITH PERIOD AND RAMP 


COMPLEX RAMP 
WIDTH .... 


47 


TABLE 


XIV. 


- 


VARIATION OF A DIFFERENTIATED VELOCITY 
SIGNAL WITH RAMP WIDTH 


52 



6 






f 






% 



% 






V 



ACKNOWLEDGEMENT 



The author would Like to thank Professor N. E. J. Boston, 
whose appreciated advice and encouragement greatly contributed 
to this thesis. The author is also grateful to Professor 
K. L. Davidson for his willing assistance. Miss S. D. Raney 
of the W. R. Church Computer Center provided much needed pro- 
gramming advice and Professor A. L. Schoenstadt provided 
mathematical assistance. Finally, the author would like to 
thank his wife, Linda, who served as part-time typist and 
proof reader. 



7 



t#’ 



<)#t 







.Its 



:t* 






i 



I . INTRODUCTION 



A. OBJECTIVE 

Atmospheric and oceanic tvirbuLence contributes greatly to 
interactions at the air/sea interface. Mathematical models 
have been developed which attempt to describe this turbulent 
boundary layer where atmosphere and ocean meet. As sensing 
systems become more sophisticated and more turbulence data 
becomes available, existing models are accordingly modified. 

An important modification was made by Kolmogorov (1962) and 
others to account for rather large fluctuations observed in 
the energy dissipation rate. Siich fluctuations were observed 
primarily in differentiated turbulence signals which emphasized 
high frequency components. Kolmogorov revised his theory in 
order to incorporate how the velocity derivative’s skewness 
(s) and coefficient of excess (c.e.) increased with the 
Reynold’s number (Re). In the energy containing region the 
parameters s and c.e. are zero since the turbulence fluctuations 
seem to have a Gaussian distribution. At small scales (high 
frequencies) they are decidedly not zero. The skewness and 
coefficient of excess of the derivative of the temperature 
fluctuations are not zero either. One might then ask, ’’Knowing 
s and c.e. from real data, is it possible to construct signals 
that resemble observed temperature signals?” Such was the, 
objective of this thesis. 
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B. GENERAL FEATURES OF TEMPERATURE SIGNALS 

The ’‘ramp” is an often observed feature in tenaperature 
fluctuations during unstable atmospheric conditions. It is 
characterized by a gradual increase in temperature followed 
by a sudden drop to an ambient level. The signal rarely falls 
below this ambient level (Figure 1). This sudden drop clearly 
distinguishes a temperature signal from other turbulence 
signals such as velocity signals whose main perturbations 
from the ambient level are usually gradual rises with 
accompanying gradual descents. It seems natural to conclude 
that these observable non-Gauss ian type features in temperature 
signals probably contribute to the measureable skewness and 
coefficient of excess. 

Electronic pattern recognition studies of atmospheric 
temperature fluctuations have been made by Sokol (1971). 

Although his study was not exhaustive, he achieved over 707o 
detection success. Another method of describing these signals 
is statistically through skewness and coefficient of excess. 

Since these statistics are zero for Gaussian noise, their 
values would represent only the non-Gauss ian character of the 
signal. Thus, it may be possible in the future to use the 
frequency and structure of the ramps as a measure of some 
meteorological conditions . 

Kaimal and Businger (1970) and Gill (1971) have demonstrated, 
over land at least, that the ramps are manifestations of 
convective plumes present near the air/earth interface. A 
schematic of the plume and its relation to the ramp, by 
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Gill (L97i), are shown in Figure 2. Taylor (1958) and Gill 
(1971) have shown that these ramps can be correlated to other 
ramps appearing simultaneously at different heights but with 
a slight offset owing to a tilt in the plume; presumably due 
to wind shear. Boston (1971) indicated that the ramps in the 
temperature signal often correlate with perturbations in 
velocity signals. Safley (1972) has related them to water 
waves and velocity signals, 

Plxjmes form a significant feature in the atmospheric 
boundary layer. Although they are present only a small 
fraction of the time they may account for the majority of heat, 
moisture, and momentum exchange. Safley (1972) found that 
the plumes accounted for 10% of the signal but 32% of the 
heat flvix. Furthermore, they cause random fluctuations in 
electromagnetic waves over line-of-s ight paths. Of recent 
interest is the effect they have on laser beams. When beamed 
over water they exhibit a random walking at the target. This 
is of considerable importance if a laser beam is to be used 
in a variety of practical applications, such as fire control 
systems . 



10 




i?fi« 








% m 

§ 



•a 

T 










11 



I 



X 



TEMPERATURE 




12 



Fig. 2. Schematic of Convective Plume According to Gill (1971). 
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II. SIGNALS 



A. ATMOSPHERIC TEMPERATURE FLUCTUATIONS 

When the atmosphere is stable the temperature fluctuations 
tend to be small because the vertical turbulent velocity com- 
ponent tends to be damped out by the vertical temperature 
gradient. When the atmosphere is neutral the fluctuations 
are also small but tend to have both positive and negative 
excursions due to the lack of restoring buoyant forces. 

Unstable conditions, on the other hand, exhibit the 
characteristic ramp features in the temperature trace (Figure 3), 
This is the condition that usually exists during the daytime 
in the lower atmosphere because the earth is heated and warms 
from below causing upward sensible heat transfer. 

B. MODELS OF ATMOSPHERIC TEMPERATURE FLUCTUATIONS 

Atmospheric turbulence for stable and neutral conditions 

is difficult to model since it is a broad band stochastic 
process. However, because of the characteristic ramp shape 
of temperature fluctuations in unstable conditions, some 
attractive possibilities for models present themselves. The 
model for the temperature signal of unstable conditions will 
consist of a basically straight line ramp signal to which 
variations will be added. 

Rather than examining spectra of the modeled atmospheric 
temperature fluctuations, statistics of the derivative of the 
model will be examined. The differentiated signal emphasizes 
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the high frequency (high wave number) region of the signal 
which is the region of most interest in the modeling process 
since measurements in the last decade clearly indicate that 
the high wave number region is decidedly non-Gaussian, Also, 
visual inspection of ramps indicates that high frequency 
fluctuations may be important in ramp structures. 
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FIGURE 3« Growth of Ramps as Atmosphere Changes from 
Stable to Unstable (Bottom to Top), 
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III. METHOD OF EXAMINATION 



A. SIGNAL 

The signal modeling was done on an IBM 360/67 digital 
computer at the W. R. Church Computer Center, Naval Postgraduate 
School, Monterey, California, The basic ingredient in each 
model is the individual ramp which will be repeated at a 
specified interval. This collection of repetitive ramps will 
comprise the model signal. Each signal will be stored in a 
10,000 X 1 vector array which represents the time axis of the 
signal. The size of the array may be varied, however, since 
the models are periodic. The distance between each array 
entry is constant and equal to the distance between samples of 
a digital recorder or converter (which is a function of the 
original sampling rate). 

The models are dimensionless although voltage is implied 
and terms such as period and width must be redefined as they 
refer to digital words. Width refers to consecutive digital 
words or array entries and is used to describe the length of 
a ramp or a ramp component. Ramp period is the number of 
consecutive digital words between the start of a ramp and the 
start of the next succeeding ramp. 

Three ramp modifications of increasing complexity will be 
examined as possible models for the temperature signal. The 
first is a "standard ramp signal" whose ramps are constructed 
of a gradual rise ending in a vertical drop (Figure 4) as 
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described by Gibson ^ aJL (L971). The variable parameters of 
this signal are (1) the period of the ramp, (2) the slope 
of the ramp, and (3) the width of the ramp. 

The second modification is an ’’asymmetrical ramp” with a 
gradual rising front slope followed by a short steep back 
slope (Figure 5), The variables in this signal are (1) ramp 
period, (2) the back slope of the ramp, (3) ramp width, and 
(4) the width of the back slope. The front slope value and 
amplitude need not be included since they are functions of 
the other variables. Specifically, 

front slope = ((back slope width) X (back slope))/ 

((ramp width) - (back slope width)), amplitude = 

(back slope) (back slope width). 

The third modification, a ’’complex ramp”, is constructed 
by inscribing nearly symmetrical consecutive triangles (teepees) 
within two envelope rays of constant but different slopes such 
that the two rays approximate a ramp function (Figure 6). 

Each succeeding slope of the inscribed teepees will be greater 
in absolute value and alternate in sign. This type of model 
is chosen because close inspection of temperature strip charts 
shows that large fluctuations of increasing amplitude is a 
commonly observed ramp feature. The variables for this 
signal are (1) ramp period, (2) slope of the upper defining 
ray, (3) slope of the lower defining ray, and (4) the width 
of the ramp. The period of the internal teepees is not listed 
as a variable because, as will be explained later, the model 
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is dimensionless and, for example, halving the period or 
doubling the frequency of the internal teepees can be accom- 
plished by simultaneously doubling the ramp period and ramp 
width. 

B. DIFFERENTIATED SIGNAL 

As previously mentioned the differentiated signal will be 
examined since it emphasizes the higher frequencies that 
contribute to the non-Gaussian nature of the signals. There- 
fore, now that the direct signal has been modeled the 
differentiated signal must be computed. The differentiated 
standard ramp signal is merely a box function whose amplitude 
is equal to the slope of the ramp (Figure 4). 

The differentiated asymmetrical ramp signal is an alternating 
asymmetrical box function. Each differentiated pulse consists 
of a broad positive box whose amplitude is equal to the front 
slope of the ramp followed by a short, highly negative box 
whose amplitude is equal to the back slope of the ramp (Figure 5). 

The signal model for the derivative of the complex ramp 

(Figure 6) is constructed by letting 

dy = us lope (1), the slope of the uppermost limiting ray and, 
dx 

dy ' = Is lope (2), the slope of the lowermost limiting ray. 
dx’ 

For representation on a digital computer the implied distance 
between successive array entries is a mesh length. Then within 
the ramps each succeeding array entry represents successive 
slopes of the inscribed teepees. However, for purposes of com- 
puting these slopes the mesh length is offset by half its 
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length and the difference between successive intersections 
of the upper slope and the lower slope is computed. From 
equations (1) and (2), 
dy = (dx)(uslope) 
dy’= (dx’)dslope). 

Let n = where n is the number of mesh lengths along the 

time axis of the signal measured from the start of the ramp. 
Hence , 

dy^ = n(uslope) 
dy’n = n(lslope). 

Therefore, assuming the first slope is positive: 
x(n) = (n+1) (uslope) - (n) (Islope) 
x(n+l)= (n+2) (Islope) - (n+1) (uslope) 
x(n+2)= (n+3) (uslope) - (n+2) (Islope) 
x(n+3)= (n+4) (Islope) - (n+3) (uslope) 
x(n+4)= (n+5) (uslope) - (n+4) (uslope) 



Note that the first and last values may have to be computed 

separately since the signal must be continuous with the ambient 

level at the beginning and end of the ramp. A sample program 

with this representation is given in a separate section. 

V 



C. STATISTICS APPLIED 

One can usually discern a temperature signal from a velocity 
signal by the ramp. However, if the spectra of these signals 
are examined in the range of about 1 Hz to 1000 Hz they are 
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remarkably similar, although their statistics such as skewness 
and coefficient of excess are quite different. As a result, 
recent turbulence literature has used these two parameters 
extensively to describe turbulence signals. Therefore, skew- 
ness and coefficient of excess will be used as modeling 
parameters. The skewness and c.e. of the models will be com- 
pared to those of actual signals. 

Skewness is a measure of the amount a signal differs from 
its mean value. The skewness of a Gaussian signal is zero. 

A distribution that is skewed to the right has a positive 
skewness and a distribution that is skewed to the left has 
a negative skewness. A symmetrical signal will have zero skew- 
ness and a signal with a high skewness value will have pulses 
that are very asjmunetrical . Coefficient of excess is a 
measure of the peaJcedness of a distribution. A distribution 
that is less peaked than normal distribution has c.e. values 
less than zero, and vice versa. The c.e. of a signal is 
highly sensitive to infrequent deviations from the mean signal 
value. Therefore, a signal with ramps spaced relatively far 
apart will have a higher c.e. value than a signal with ramps 
spaced close together. These higher order statistics emphasize 
the higher frequencies v;hich are non-Gaussian. However, since 
they are exponential functions , they require a longer record 
to achieve stable results. Therefore, the use of statistics 
of orders higher than c.e. is somewhat limited. 

Skewness (s) is defined as : 

3 9 3 /2 

s = e /(e"^) , where e is the variation of the 

signal from its mean value. 
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Coefficient of excess (c.e.) is defined as: 

_ 4 4 9 2 

Often one sees e /(e ) , the ratio 

(e^ 

of the fourth moment to the second moment, used as a statistic 
and referred to as "kurtosis," but c.e. will be used here 
because it is zero for a Gaussian signal. 

The skewness and c.e. for differentiated temperature 
signals have been computed by several people (Table I). Al- 
though the one thing the statistics seem to have in common is 
their variability, representative figvires can be chosen. For 
temperature signal a skewness value between 0.2 and 0.5 seems 
common and likewise a value for c.e. is about 30. Compare 
these values to representative velocity signal values of -0.5 
and 15 for s and c.e. respectively. 
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DIFFERENTIATED ASYMMETRICAL RAMP SIGNAL 
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UPPER DEFINING 
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FIGURE 



TABLE I 



SKEWNESS AND COEFFICIENT OF EXCESS OF 
DIFFERENTIATED TEMPERATURE SIGNALS 



INVESTIGATOR 


SKEWNESS 


C.E. 


BOSTON (1971) 


-0.013 


20.3 




0.195 


17.1 




0.140 


21.4 




0.002 


22.6 




-0.082 


13.6 




0.596 


22.8 




0.203 


12.2 




0.307 


13.4 




0.315 


32.0 




0.171 


13.9 




0.673 


30.7 




0.039 


15.5 




0.286 


18.1 




0.251 


11.0 




0.257 


12.3 




0.051 


6.0 


GIBSON al. (1970) 


0.40 


40 




0.72 


35 




0.48 


34 




0.72 


22 


GURVICH (1967) 




15 - 
14001 
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IV . RESULTS 



The standard ramp signal proved to be a poor representation 
of the temperature signal. For an appropriate c.e, the skew- 
ness was positive but much too large; about five. This agrees 
with the results obtained by Gibson e_t.al. (1971) from a 
similar model. The asymmetrical ramp was no improvement. With 
this signal a large positive c.e. was accompanied by a large 
negative skewness. This signal proved to be of some value 
apart from its intended use and is discussed in Section VI. 

The complex ramp signal showed more favorable results. Not 
only can compatible skewness and c.e. be found but the signal 
can be adjusted so that skewness and c.e. vary almost indepen- 
dently. A more detailed analysis of the individual signals 
follows. It should be kept in mind that while a direct signal 
model is the final product, the different model theories were 
tested by the statistics of the differentiated signal. There- 
fore, the following results apply to the statistics of the 
differentiated signal. 
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TABLE II 



VARIATION OF DIFFERENTIATED RAMP SIGNAL 
WITH PERIOD (FIGURE 7) 



PERIOD 


RAMP SLOPE RAMP WIDTH 


S. 


C.E. 


10 




0.408 


-1.83 


20 




1.45 


0.250 


30 




2.15 


2.65 


40 


2.0 4.0 


2.63 


4.90 


80 




3.98 


13.8 


100 




4.70 


20.0 


200 




6.48 


40.0 


400 




8.49 


70.0 


Note 


here that skewness and c.e. both 


increase 


with 



increasing pericxi but they never simultaneously achieve the 
correct values required by the models. 
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FIGURE 7. Variation of Standard Ramp Signal with Period 
(Ramp Slope = 2,0, Ramp Width = 4.0) 



TABLE III 



VARIATION OF DIFFERENTIATED RAMP SIGNAL 
WITH RAMP SLOPE (FIGURE 9) 



PERIOD 


RAMP SLOPE 


RAMP WIDTH 


S.' 


C.E. 




0.1 




2.67 


5.12 




0.5 




2,67 


5.12 


40 


1.0 


4.0 


2.67 


5.12 




5.0 




2.67 


5.12 



S. and c,e. arc independent of ramp slope. 
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TABLE IV 



VARIATION OF DIFFERENTIATED RAMP SIGNAL 
WITH RAMP WIDTH (FIGURE 8) 



PERIOD 


RAMP SLOPE 


RAMP WIDTH 


S. 


C.E. 






2.0 


4.13 


15.1 






3.0 


3.23 


8.41 


40 


2.0 


4.0 


2.67 


5.11 






5.0 


2.27 


3.14 






6.0 


1.96 


1.85 






7.0 


1.71 


0.929 



Here the results are similar to those obtained with a 
varying period. Although values of s and c.e. can be achieved 
with the standard ramp signal that ai'e respectively commensurate 
with measured values from actual signals, the values of s and 
c.e, are not simultaneously compatible with measured signals 
and hence this signal is discarded as an appropriate model. 
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FIGURE 8. Variation of Standard Ramp Signal 
with Ramp Width (Period = 40, Ramp 
Slope = 2.0) 
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FIGURE 9. Variation of Standard Ramp Signal 
with Ramp S-lope (Period = 40, Ramp 
Width = 4.0) 



31 



% 






TABLE V 



VARIATION OF DIFFERENTIATED ASYMMETRICAL 
RAMP SIGNAL WITH PERIOD (FIGURE 10) 



PERIOD 


BACK 

SLOPE 


FRONT 

SLOPE 


RAMP 

WIDTH 


BACK 

SLOPE WIDTH 


S. 


C.E, 


10 










-3.77 


13.2 


30 










-4.64 


21.6 


50 










-6.00 


38.0 


100 


-1.5 


0.167 


10 


1 


-8.48 


78.9 


300 










-14.7 


243 


500 










-18.9 


403 



Although c.c. values are achieved that agree with actual 
signal values, the skewness values are not only too large in 
absolute value but negative instead of the small positive 
values desired. 
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SKEWNESS AND COEFFICIENT OF EXCESS 




FIGURE LO. Variation of Asymmetrical Ramp Signal 
with Period 

(Back Slope = -1.5, Ramp Width = 10, 
Back Slope Width = 1) 
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TABLE VI 



VARIATION OF DIFFERENTIATED ASYMMETRICAL RAMP 
SIGNAL WITH BACK SLOPE (FIGURE 12) 



PERIOD 


BACK 

SLOPE 


FRONT 

SLOPE 


RAMP 

WIDTH 


BACK SLOPE 
WIDTH 


S. 


C.E 




-0.10 


0.011 






-8.43 


78.1 




-0.50 


0.056 






-8.43 


78.1 


100 


-1.0 


0.10 


10 


1 


-8.43 


78.1 




-1.5 


0.144 






-8.43 


78.1 




-1.9 


0.189 






-8.43 


78.1 



Once again the slope of the signal has no effect on the 
modeling statistics. 
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* 



VARIATION OF DIFFERENTIATED ASYMMETRICAL RAMP SIGNAL 
WITH RAMP WIDTH (FIGURE LI) 



PERIOD 


BACK 

SLOPE 


FRONT 

SLOPE 


RAMP 

WIDTH 


BACK 

SLOPE WIDTH 


S. 


C.E. 






0.30 


6 




-7.30 


67.0 






0.25 


7 




-7.72 


70.8 






0.21 


8 




-8.02 


73.8 


100 


-1.5 


0.19 


9 


1 


-8.25 


76.2 






0.17 


10 




-8.43 


78.1 






0.08 


20 




-9.23 


87.3 






0.05 


30 




-9.49 


90.4 



These statistics are not compatible with measured values. 
The skewness is too large in absolute value and negative in 
sign. 
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SKEVVNESS AND COEFFICIENT OF EXCESS SKEWNESS COEFFICIENT OF EXCESS 




FIGURE 11: VARIATION OF ASYMl'IETRICAL 
(PERIOD “ 100, BACK SLOPE 



RAMP SIGNAL ’WITH RAMP WIDTH 
=-1.5, BACK SLOPE V/IDTH = 1) 




FIGURE 12: VARIATION OF ASYMMETRICAL RA^IP SIGNAL WITH BACK SLOPE 
(PERIOD » 100, RAl-D? vm)TH “ 10, BACK SLOPE WIDTH = 1) 
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TABLE VIIIo 



VARIATION OF DIFFERENTIATED ASYMMETRICAL RAMP 
SIGNAL WITH BACK SLOPE WIDTH (FIGURE 13) 



PERIOD 


BACK 

SLOPE 


FRONT 

SLOPE 


RAMP 

WIDTH 


BACK SLOPE 
WIDTH 


S. 


C.E. 






0,167 




1 


-8.48 


78.9 






0.375 




2 


-4.77 


29.8 


100 


-1,5 


0.643 


10 


3 


-2.77 


14.8 






1.0 




4 


-1.30 


S.7i 






1.5 




5 


0.0 


7.0 



Note that the skewness finally becomes non-negative (zero) 
but the c«,e. is too low and the shape of the signal is that of 
a teepee, not a ramp. Therefore, due to the negative skewness 
inherent in this signal, it is discarded as a possible model. 
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VARIATION OF ASYMMETRICAL RAMP SIGNAL VTETH BACK SLOPE \VIDTH 
(PERIOD « 100, BACK SLOPE » -1,5, RAMP ^vTDTH = 10) 

FIGURE 13. 
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TABLE IX 



VARIATION OF DIFFERENTIATED COMPLEX RAMP SIGNAL 
WITH PERIOD (FIGURE 14) 



PERIOD 


UPPER BASIC 
SLOPE 


LOWER BASIC 
SLOPE 


RAMP S. 

WIDTH 


C.E, 


40 






0.046 


-1.27 


80 






0.064 


0.45! 


100 






0.072 


1.32 


200 






0.102 


5.64 


300 


0.8 


-0.4 


40 0.124 


9.71 


400 






0.145 


14.3 


600 






0.175 


22.4 


800 






0.201 


30.2 


1000 






0.229 


40.2 



This signal closely approximates actual measured values 
for periods in the range of 800 to 1000. This appears to be 
an excellent model for typical atmospheric temperature signals 
during unstable conditions . Now the effect of the remaining 
parameters will be examined. 
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200 



Uoo 600 

PERIOD 



FIGURE L4 : VARIATION OF C0^^PLEX RAMP SIGNAL WITH PERIOD 
(UPPER BASIC SLOPE - 0.8, LO\^ER BASIC SLOPE = -O.U, RAI'IP WIDTH » lO) 
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TABLE X 



VARIATION OF DIFFERENTIATED COMPLEX RAMP SIGNAL 
WITH UPPER BASIC SLOPE (FIGURE 15) 



PERIOD UPPER BASIC 

SLOPE 


LOWER BASIC 
SLOPE 


RAMP 

WIDTH 


S. 


C.E. 


0.4 






0.208 


30.6 


0.8 






0.201 


30.2 


1.2 






0.175 


30.2 


800 1.6 


-0.4 


40 


0.153 


30.2 


o 

• 

CM 






0.134 


30.2 


2.4 






0.119 


30.2 


50 






0.008 


30.7 


The coefficient of 


excess remains 


almost 


constant 


while 



the skewness ranges over values not unlike those encountered 
in real signals « 
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FIGURE 15 :VAKIATION OF COMPLEX RAMP SIGNAL V/ITH UPPER BASIC SLOPE 

(PERIOD * 800, LOVffiR BASIC SLOPE = -O.U, HAI-IP VvTDTH = Uo) 



42 



TABLE XI 



VARIATION OF DIFFERENTIATED COMPLEX RAMP SIGNAL 
WITH LOWER BASIC SLOPE (FIGURE 16) 



PERIOD 


UPPER BASIC 
SLOPE 


LOWER BASIC 
SLOPE 


RAMP 

WIDTH 


S. 


C.E. 






0.6 




-10.3 


227 






0.4 




- 2.29 


58.2 






0.2 




- 0.488 


34.6 


800 


0.8 


0.0 


40 


0.0 


30.7 






-0.2 




0.153 


30.2 






-0.4 




0.201 


30.2 






-0.6 




0.211 


30.4 


The 


values of s and 


c.e. obtained 


when the 


lower basic 



slope is positive fluctuate greatly but when the lower basic 
slope becomes negative the c.e. becomes almost constant, about 
30, and the skewness stabilizes near 0,2. Therefore, it is 
concluded that the model of a temperature signal has a lower 
basic slope less than zero. 
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16: VARIATION OF COMPLEX RAMP SIGNAL WITH LOWER BASIC 

SLOPE (PERIOD = 800, UPPER BASIC SLOPE = 0.8, 

RAMP WIDTH = 40) 
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TABLE XII. 



VARIATION OF DIFFERENTIATED COMPLEX RAMP SIGNAL 
WITH RAMP WIDTH (FIGURE 17) 



PERIOD 


UPPER BASIC 
SLOPE 


LOWER BASIC 
SLOPE 


RAMP 

WIDTH 


S. 


C.E. 








40 


0.201 


30.2 








50 


0.145 


23.8 








60 


0.111 


19.5 


800 


0.8 


-0.4 


70 


0.089 


16.3 








80 


0.073 


14.0 








90 


0.061 


12.1 








100 


0.052 


10.6 



The effect of varying the ramp width has much the same 
effect as varying the period. Decreasing the ramp width 
causes s and c.e. to increase; likewise an increase in period 
has similar results. These two adjustments however do not 
produce equal changes in s and c.e. For example, in TABLE IX 
there is a signal with skewness of 0.145 and a c.e. of 14,3, 
while from TABLE XII a signal with skewness of 0.145 has a 
c,e, of 23.8. 
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FIGURE 17: VARIATION OF COMPLEX RAMP SIGNAL VJITH R/iMP WIDTH 

(PERIOD “ 800, UPPER BASIC SLOPE - 0.8, LOVJER Bi\SIC SLOPE - -0.4) 
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TABLE XIII. 



VARIATION OF DIFFERENTIATED COMPLEX RAMP SIGNAL 
WITH PERIOD AND RAMP WIDTH 



PERIOD 


UPPER BASIC 
SLOPE 


LOWER BASIC 
SLOPE 


RAMP 

WIDTH 


S. 


C.E. 


200 






10 


0.702 


28.1 


400 






20 


0.389 


30.3 


800 


0.8 


-0.4 


40 


0.201 


30.2 


1600 






80 


0.099 


28.5 


3200 






160 


0.047 


24.8 



TABLE XIII is Listed for the information of the reader and 
as an aid to those who might want to model temperature signals. 
In each successive entry in TABLE XIII the period and ramp 
width have been doubled. As previously noted, this has the 
effect of doubling the frequency of the internal teepees. 

Also, it is a means of changing the skewness while keeping 
the c.e. relatively constant. 

The preceding tables have not been an attempt to list all 
the variations possible with the different models, but to 
examine how various parameters affect the models and to 
determine which model, if any, is the best representation of 
atmospheric temperature signals. As previously noted, the 
complex ramp signal was selected as the best representation. 
Obviously there exists an infinite number of possible variations 
and given the actual values of a signal together with the 
previous tables one should be able to model most of the 
atmospheric temperature signals encountered in an unstable 
atmosphere . 
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The complex ramp provides a coarse and a fine "tuning'* 
mechanism for the model. The ramp width and/or the period 
can be adjusted to the desired c.e. value with an s value 
near that desired and then the skewness tuned by adjusting 
the upper basic slope while the coefficient of excess remains 
constant. Additionally, the model of the differentiated 
complex ramp (Figure 6) is remarkably similar to the actual 
differentiated signal (Figure 1). 
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V. CONCLUSIONS 



The results of this study indicate successful accoiaplish- 
ment of the objective; which was to model a temperature signal 
from its statistical parameters, skewness, and coefficient 
of excess. Furthermore, additional insight was gained into 
the actual structure of temperature signals. 

Although the study was not exhaustive it does indicate 
that one cannot attribute the observed values of skewness and 
coefficient of excess in the differentiated temperature signal 
entirely to the external dimensions and frequency of the ramps. 
Instead one must look at the internal structure of the ramp 
itself. The ramp is composed of a great many fluctuating 
high frequency saw-teeth which gradually increase in amplitude 
giving an overall impression of a relatively gradual rising 
ramp. The fluctuations are presumably due to mixing or 
entrainment by the diffuse leading edge of the manifesting 
plume. The saw- teeth actually extend below the ambient 
temperature level. This is not readily apparent when examining 
temperature strip charts but can be attributed to the mechanical 
response of the recorder and/or sensing system. The internal 
saw-teeth cause the variation in skewness while the width or 
frequency of the entire ramp causes variation in c.e. 

Therefore, given the skewness and coefficient of excess of a 
temperature signal, a model can be constructed with the exact 

statistics of the original signal. 

! 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 



A natural extension of this thesis would be an attempt 
to model actual signals. This would probably require not 
only an examination of strip chart recordings but also an 
investigation of the signal's digital representation. 

While the asymmetrical ramp signal proved unsatisfactory 
as a temperature model, in an extreme variation it's 
statistics were remarkably similar to those of a velocity 
signal. When the back slope width was only slightly smaller 
than half the entire ramp width a small negative skewness 
and a c.e. of near 20 was noted. Several of these results 
are listed in TABLE XIII. Except for a shape that more 
resembles a teepee than a ramp, this signal should behave 
exactly like the previously noted asymmetrical ramp signal. 

An inspection of a velocity signal indicates that this teepee- 
like representation is not a bad resemblance. However, it 
should also be noted that the same results can be achieved 
with a complex ramp signal whose lower basic slope is greater 
than zero but necessarily less than the upper basic slope 
(TABLE XI). Perhaps a more accurate model would be a combi- 
nation of the two, where high frequency fluctuations are 
inscribed between two asymmetrical ramps similar to those 
noted above such that the ramps have the same widths but 
different heights, i.e., front slopes and back slopes. 
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If this model proved successful it would hopefully 
inherit the characteristic of the complex ramp whereby the 
skewness and c.e. could be adjusted independently. 
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TABLE XIV 



VARIATION OF A DIFFERENTIATED VELOCITY 
SIGNAL WITH RAMP WIDTH 



PERIOD 


BACK 


FRONT 


RAMP 


BACK SLOPE 


s. 


C.E. 




SLOPE 


SLOPE 


WIDTH 


WIDTH 










L.5 


4 




0.0 


22.0 






L.O 


5 




-L.53 


20.3 


LOO 


-1.5 


0.75 


6 


o 

• 

CM 


00 

• 

CM 

1 


22.0 






0.60 


7 




-3.59 


24. L 
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THIS PROGRAM COMPUTES TH= SKSWMHSS (SKEW) AND CQHFFICISMT 
OF 'XCcSS (C~) FOP. A COMPLEX R4MP SIGNAL WITH A PERIOD OF 
8GC, AN UPPER BASIC SLOPE OF 0.8, A LCW2R BASIC SLOPE OF 
-0.4 AND A RaMP WIDTH OF 40. PARAMETERS BOXED WITH • S APS 
THE VARIABLES OF TH<^ SIGNAL. A DO LOOP INS.RTSO BEFORE 
STATEMENT 5 AND LOOKED THROUGH STATEMENT 90 WILL ALLOW ONE 
OR MORE OF THE VARIABLES TO BE INCREMENTED. 

INITIALIZE DATA 

REAL*4 LSLOP? 

DIMENSION X(IOOCO) 

5 DO 10 1 = 1, 10000 
10 X(I)=0.0 
N0=10000 
M = 1 

CODE = lC*^nO 

*' I NT =80 C 
^ USLOPS=O.S 
•? LSL0PE=-0.4 * 

* VIIDTH=40 

IND=0 

Nl=2 

N2=WIDTH/2 

COMPUTE COMPLEX SIGN.AL 

DO 50 K=1,1COOO 

CCMDIJTS POSITIVE SLOPES 

X( 1 + IND ) = 2.-*USL0PE 
CO 2 5 J = N1 ,N2 

25 X( 2*- J-l + IND) =2.^ J-USL0P*^-(2.^- J-1 litLSLOPE 
CCMPUTE NEGATIVE SLOPES 
N3=N2-1 

X(2 + 1ND)=3.-'LSL0PS-2.:*:USL0PE 
DO 35 J=N1 ,N3 

35 X( 2«J + IND) = ( 2.'"J+1) •LSLOPE-2.^- JAUSLOPE 
X (2*M2 + IN0 ) = -2 .>^N2 ■■USLOPE 
INO=IMO+IMT 

CH^CK IF SIGNAL ARRAY IS FILLED 

IF( 2‘N2+IND.GT. 10000) GO TO 75 
50 CONTINUE 

75 CALL MI SR( NO ,M,X ,CODE, XRAR,STC, SKfW,CE,R ,N,A ,8.S , lER ) 
80 WRIT^ (6,20C) !NT,USLnP",LSLOPE, WIDTH, SKEW, C5 
85 CONTINUE 
90 WRITE (6,300) X 

200 FORMAT ( • C • , ' INT = • , I 5 , 2X , 'USL CP H= • , F7 . 3, 2X, *L SL OP^ = • , 
1F7. 3,2X,'WIDTH = ' ,F5.1,2X,' SK?W=' , F8 . 3 , 2X , ' CE = ' , =8 .3 , / ) 
300 FORMAT (IX, 'DATA SIGNAL FOLLOWS './, 1C0( lOF 10. 3/ ) ) 

STOP 

END 
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SUBROUTINE MI SR 
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